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We consider here the non-adiabatic energy transfer dynamics for a model bi-chromophore system consisting of a perylenedi- 
imide unit linked to a ladder-type poly-(para-phenylene) oligomer. Starting from a semi-empirical parameterization of a model 
electron/phonon Hamiltonian, we compute the golden-rule rate for energy transfer from the LPPP5 donor to the PDI acceptor. 
^ Our results indicate that the non-adiabatic transfer is promoted by the out-of-plane wagging modes of the C-H bonds even though 
f— s theses modes give little or no contribution to the Franck Condon factors in this system. We also predict a kinetic isotope effect 
{S( ofk^/k^ = 1.7 — 2.5 depending upon the temperature. 
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1 Introduction 

Electronic energy transfer between donor and acceptor units 
provides the basic energy transport mechanism for optical- 
electronic devices and photo synthetic systems in nature. For 
the case of separated donor/acceptor species, one typically as- 
sumes that the off-diagonal coupling between states, /, is inde- 
pendent of the internal vibrational motions of the two species 
and that the internal motions of each species are independent 
of each other allowing one to treat them as "separate baths." 
This ansatz is useful since it allows one to compute transfer 
rates based upon the spectral overlap between isolated donor 
j and acceptor states. This, along with a other assumptions such 
as that the coupling can be estimated by using the transition 
dipole moments allows one to write the transfer rate as 

k DA = ^2 / ^a(co)/ d (co)Jco (1) 

27T/T JO 

where Id (to) is the fluorescence spectrum of the donor and 
£a(co) is the (normalized) absorption spectrum of the accep- 
tor. The coupling / is the dipole-dipole coupling between tran- 
sition moments which scales with as l/R 6 in the separation 
between donor and acceptor species. This is valid only when 
R is large compared to size of the chromophores themselves. 
With in the Forster model, energy transfer is becomes efficient 
when there is sufficient overlap between the emission spec- 
trum of the donor and the absorption spectrum of the acceptor. 
However, in a number of interesting cases energy transfer 
can be very efficient even in the absence of significant spec- 
tral overlap™. This scenario can arise in cases where the 
surrounding media appears to participate in the energy trans- 
fer process by retaining memory of the donor state long after 
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LPPP5 



Fig. 1 Chemical structure of donor- acceptor dyad: LPPP5 
(ladder- type poly-(/?-phenylene)) PDI (perylenediimide). The 
highlighted bonds (in red) define the dihedral angle, (]), between the 
donor and acceptor moieties. For the ground, donor, and acceptor 
states § gs = 59.6°, § D = 44.2°, and § A = 59.1°, respectively. 



the transfer has occurred. Here, it is often useful to invoke 
within the model a common "bath" or "shared" sets of modes 
between the donor and acceptor species such that the common 
mode not necessarily appearing in the absorption or emission 
spectra participates in the energy transfer event 5 . However, 
there remains some ambiguity in describing this mode and its 
coupling to the electronic transition. 

Our goal in this paper is to devise a systematic approach 
for determining the types of modes that participate in phonon- 
assisted energy transfer process and to test this approach on 
a series of donor-acceptor species known to exhibit efficient 
energy transfer. Our starting point will be from a quantum 
chemical perspective which will allow us to parameterize our 
theoretical description based upon molecular considerations 
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and systematically improve our description using the appro- 
priate levels of quantum chemical theory. 

In this study, we focus on supramolecular donor- acceptor 
(D-A) unit 6 , which is made up of a donor molecule 
"LPPP5" (Ladder-type poly-(para-phenylene)) and an accep- 
tor molecule "PDF (perylenediimide) linked by a single cova- 
lent bond. The chemical structure of the D-A unit are shown 
in Fig.[T] Although donor and acceptor are connected by a sin- 
gle covalent bond, quantum chemical calculations have shown 
that the molecular orientation of LPPPI and PDI are tilted 
along the single bond approximately by 60° with respect to 
each other. LPPP and PDI has been extensively studied and 
has wide application potential 7 . LPPP is a p-type electron 
donating semiconducting polymer which is subject to a par- 
ticularly low degree of disorder, as the full ladder structure 
along the backbone blocks chain coiling and bending enforc- 
ing a rigid-rod type of geometry 7 8 . A rigid coplanar structure 
of conjugated ladder- type polymers are very suitable for light- 
emitting applications which enhances the conjugation, carrier 
mobility, and luminescence intensity 7 . PDI polymers show 
an n-type, electron conducting behavior and serve as electron- 
acceptor material 9 10 . Covalently linked molecular D-A units 
can be seen as model compounds for D-A blends used in op- 
toelectronic devices. After photoexcitation of D-A unit, exci- 
tation energy transfer occurs from an excited donor (D*A) to 
an acceptor molecule (DA*) to a degree determined by dipole- 
dipole interactionLQ-l 

2 Theoretical approach 

2.1 Non-adiabatic Hamiltonian 

For the co-joined species, we can not make a clear separation 
between donor and acceptor species since there is the possi- 
bility of strong electronic mixing and sharing of vibrational 
modes. Our approach here will be to work within the elec- 
tronic eigenstate basis and assume that we can treat the elec- 
tron/phonon coupling as linear in the vibrational coordinates, 
so that the entire molecule can be treated within a linearized 
Born-Oppenheimer approximation. We also assume that the 
normal modes of the ground state are a good description of 
the modes in the excited states. To begin, we shall expand 
the electronic Hamiltonian about the ground state equilibrium 
geometry of the molecule. 



H(q)=H(0)+VHq+^V 2 Hq 



(2) 



Where vector q denotes all the normal mode coordinates of 
the supramolecular system 6 at the ground state equilibrium 
geometry (q = 0). It represents mass-scaled normal coordi- 
nate local to the entire donor-acceptor unit with vibrational 
frequency CO. We compute electronic eigenstates along the q 



while satisfying 

H(q)\y a (q))=E a (q)\\\f a (q)). 

Thus, the matrix elements of H(q) can written in a primitive 
adiabatic basis as^L^ 



H(q) 



E a (0)+ga-q VqHab-q 

V q H ba -q E b (0)+g b -q 



1 



coV 



(3) 



Where, E a j,(Q) are the vertical energies at the ground state 
equilibrium geometry (q = 0) . g a is the force directed along 
the vector from the ground state equilibrium geometry to adi- 
abatic minima of excited state \\f a (q) . Similarly, g b is the force 
directed along the vector from the ground state equilibrium 
geometry to the adiabatic minima of excited state \\f b (q) . We 
can interpret g a and g b as being two non-orthogonal vectors 
originating at the ground state equilibrium geometry of the 
supramolecular system and pointing towards the equilibrium 
geometry of the respective excited states 14 . 

Within a harmonic approximation, we can determine the di- 
agonal terms by taking the energy derivative of E a along a 
given normal mode direction. 



A (Yfl |H| Vfl )=^ 

aq n aq n 



(4) 



qn=0 



The on-diagonal couplings can be determined numerically 
from quantum chemistry as the energy gradient of an excited 
state energy taken at the ground- state equilibrium geometry. 
Thus, g aan are derived from the reorganization energy and shift 
in the geometry from ground state (q = 0) to the excited state 
(a or b). In the harmonic approximation, the equilibrium ge- 
ometry of the molecule in a given electronic excited state is 
shifted along each normal mode relative to the ground state 
equilibrium geometry by g a an/®n with energy 



aan 

2 ' 



CO: 



The off-diagonal terms, VH ab ■ q are the non-adiabatic cou- 
plings between electronic eigenstates \\f a (q) and \\f b (q) at the 
ground state geometry (q = 0) belonging to the adiabatic po- 
tentials at E a (q) and E b (q). Since \\f a (q) and y b (q) are elec- 
tronic eigenstates of H(q) for all q, V q H ab (q) = 0. Conse- 
quently, close to the ground- state geometry, the non-adiabatic 
coupling can be rewritten as 



V q H ab -q = (y a (0)\V q H(0)\y b (0)).q 
= (E a -E b )(y a \V q \y b ).q. 



(5) 



This is often referred to as the "off-diagonal" Hellmann- 
Feynman theorem 15 . We next replace the gradient operator 
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Table 1 Quantum chemical (AM1/CI(S)) data for LPPP5-PDI dyad. E n is the vertical excitation energy and bE n is the reorganization energy 
of each state. r^A is the length of the connecting bond and § is the dihedral angle between LPPP5 and PDi moieties for each optimized 
electronic state. 

state designation E n (eV) dE n (meV) ?da(A) § 



So 


DA 





- 


1.38 


59.6° 


Si 


DA* 


2.68 eV 


82.7 


1.40 


44.2° 


s 2 


D*A 


2.98 eV 


142.6 


1.38 


59.1° 



S 2 -► Si transition moment (Debye): ^12=0.063 x + 0.078 y + 0.007 z 



with the phonon momentum operator, p = —hiV q . Then, we 
recognize that if we write the molecular Hamiltonian in the 
mass- scaled coordinates as 



H-mol ' 



1 



-H(q) 



(6) 



then [H moh q] = 
write 

VqHab ' q - 



ihp. Thus, V q 



[H mo i,q]/h 2 and we can 



(E a -E b ) 



{Va\[H mo uq]Wb)-q. 



(7) 



Next, we assume that (v|/ fl |/? 2 |\|/^) 
that 



and can be ignored so 



VqHab ' q : 



(Eg-E h f 

h 2 e 



(Pab • q) 



(8) 



where p ab = e(\\f a \q\\\f b ) is the electronic transition dipole mo- 
ment between states \\f a and \\f b computed at q = 0. Here we 
are using the normal mode coordinates as a general basis for 
the position operator that can act on on the electronic degrees 
of freedom. A simple justification for this is that the Hamil- 
tonian H must be in the totally symmetric irreducible repre- 
sentation, thus the components of VH must be in the same 
irreducible representation as x, y, and z respectively. Thus, 
the non-adiabatic coupling can be approximated by taking the 
projection of the electronic coupling between states \\f a and \\f b 
and projecting this along the displacement vectors for normal 
mode q. Rewriting g abn = VH ab ■ q and p abn = p ab ■ q, which 
gives 



§abn 



(Eq-E b ) 2 

h 2 e 



H-abn 



(9) 



as the non-adiabatic (off-diagonal) coupling for n th normal 
mode. Hence, the g abn is calculated by projecting the dipole- 
transition moments p ab between excited states onto the mass 
weighted normal mode vector q. It is important to recognize 
that within the molecular/non-adiabatic model, nuclear mo- 
tions that lead to the geometric distortions of the molecule in a 



given electronic excited state may not necessarily be the same 
set of modes that couple the two electronic states. 

Let us now write H in terms of phonon operators q n = 
^^/^(^(al -\- a n ) and define G a bn — ^/h/^ngabn- Here, all 
the G a bn parameters are in units of energy rather than units of 
force. 



H = Y, E a\ a )( a \ + J^GabnW)(b\(al-\-an) 



abn 



Y,n® n (aia n + -) 



(10) 



Here |a)'s denote electronic states with vertical energies e a , 
a\ and a n are the creation and annihilation operators for the 
normal mode n with frequency co n , and G abn are the coupling 
parameters of the electron-phonon interaction which we take 
to be linear in the phonon normal mode displacement coordi- 
nate. 

We can separate H into a part that is diagonal with respect 
to the electronic degrees of freedom, 

H = J^E a \a)(a\+Y,GaanW)(a\(al + a n ) 

a an 

+ Y< hay n(aia n +-) (11) 



and an off-diagonal part V 

V = Y,Gabn\a)(b\(al + a n ), 



(12) 



abn 



where the prime at the summation sign indicates that the terms 
with a = b are excluded. This separation is useful for the fol- 
lowing two reasons. First, in many systems only off-diagonal 
coefficients G a u are small compared to G aan - Hence, V can be 
treated as a perturbation. Second, for many cases of interest, 
the initial density matrix commutes with Hq. In this case, the 
separation gives simpler forms of the master equations. 



2.2 Obtaining 
chemisty 



model parameters from quantum 



One of the advantages of our approach is that one can in 
principle arrive at a complete parameterization of our model 
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S2 transition density charges 



Si transition density charges 



Fig. 2 Molecular orbitals and energies of the co-joined LPPP5-PDI 
system along with the transition charges for the S — > S\ and 
S — > $2 transitions. 



Hamiltonian from quantum chemical considerations. For the 
LPPP5-PDI di-chromophore unit shown in Fig. [T] we used 
the semi-empirical AM1/INDO model as implemented in the 
AMPAC package 1617 . This approach is robust for system 
such as this and gives reliable energetics and geometry with 
a reasonable amount of computational overhead. We first 
performed ground- state optimization and normal mode anal- 
ysis and then performed configuration interaction (CI) calcu- 
lations to obtain the optimized geometries of the two lowest 
singlet electronic excited states. We also obtained the transi- 
tion dipole moment between excited states using this proce- 
dure. Relevant data from these calculations are presented in 
Table[Q 



« B k 
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(o (cm -1 ) 
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(a) 



(b) 



Fig. 3 (a) Contribution from each normal mode to the 
reorganization energy of the donor (black: solid line) and acceptor 
(dotted line) states for the LPPP5-PDI di-chromophore system, (b) 
Magnitude of transition dipole moment vector ~pL a \ } as projected on to 
each normal mode between states S\ and S 2 




Fig. 4 Projection of DA transition moment, JIdAj onto atomic 
displacements for PDI-LPPP5 dyad. 



For LPPP5-PDI dyad, the S 2 donor state lies 2.98 eV above 
the ground state and is largely localized on the LPPP5 end of 
the molecule. The S\ (acceptor) state is lower at 2.68 eV. Both 
are optically coupled to the ground states. If the system were 
forced to be planar, the ^-conjugation would extend across 
the entire molecule. The two components of this molecule 
are linked by a bond connecting C atoms participating in the 
7T conjugation of both moieties. In the ground and relaxed 
S 2 states, the dihedral angle between conjugated domains is 
nearly 60°. In the S\ (acceptor) state the C-C bond connect- 
ing the two moieties increases slightly to 1 .40 A and the sys- 
tem becomes slightly more planar ((|) = 44.2°). In the elec- 
tronic ground state, the dihedral angle formed by the respec- 
tive molecular planes is close to 60° suggesting that the 7T sys- 
tem should be localized on the donor and acceptor moieties. 
In Fig. [2] we show the HOMO-1 through LUMO+1 molecular 
orbitals for the dyad along with the associated transition den- 
sities for the S — » S n transitions. By and large, the HOMO 
and LUMO orbitals are localized on the PDI side of the dyad, 
although there is a significant amount of 7T amplitude leak- 
ing though the linkage over to the LPPP5 end. Likewise, the 
LUMO-1 and HOMO+1 are mostly localized on the LPPP5 
end, but one can clearly see significant leakage over to the 
PDI end. 

At the top of Fig. [2] we show the corresponding transition 
densities for the S -> Si (2.68 eV) and S -> S 2 (2.98 eV) 
transitions. Since the AMI model is based upon the zero- 
differential overlap approximation, the transition densities are 
represented as atom centered charges rather than spatial densi- 
ties. The transition moment is obtained by multiplying the lo- 
cal charge by the atomic coordinate vector and summing over 
all atoms. By and large, the vertical transitions are localized to 
the respective donor and acceptor sides of the molecule with 
some "leakage" across the covalent linking bond to the other 
side of the molecule. This weak ^-communication could in 
part account for the efficient non-adiabatic relaxation in this 
system. 

In Fig. [3^ we show the contribution to the total adiabatic 
reorganization energy from each normal mode for both the 
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donor and acceptor states. As discussed above, these are pro- 
portional to the magnitude of the diagonal contribution to the 
electron/phonon coupling. For the S\ state, the largest contri- 
bution to the reorganization energy comes from a set of low 
frequency modes (c.a. 30 cm -1 ). These correspond to the 
frustrated dihedral torsional motion at the linking point be- 
tween the two parts of the molecule. This feature is not present 
in the coupling spectrum for the other state. Common to both, 
however, are contributions from modes around 1600 cm -1 and 
1800 cm -1 . These correspond to the C=C bond stretching 
modes and are consistent with the the vibronic fine- structure 
peaks observed in most conjugated polymer systems^^TCS 

Fig. [3]) shows the projection of the transition moment be- 
tween the S\ and S2 states onto the normal modes of the 
system. As noted above, these are related contributions to 
the non-adiabatic couplings between the two states. Here we 
note two distinct contributions to the non-adiabatic coupling. 
These correspond to the wagging (at 1200 cm -1 ) and stretch- 
ing (at 3100 cm -1 ) of C-H bonds attached to the conjugated 
rings of the system. In Fig. [4] we illustrate this graphically 
by drawing the projection of the transition moment on to the 
individual atomic displacements of the system. The relative 
length of each vector indicates the relative component of par- 
ticular atomic displacement along the total transition moment 
when summed over normal modes. Surprisingly, the modes 
that contribute strongly to the adiabatic reorganization of each 
state give little contribution to the non-adiabatic coupling be- 
tween the two states. This is surprising since one expects that 
the states with the largest electron-phonon coupling would 
give both the largest contribution to both the reorganization 
and the state- to- state transitions. 




0.005 0.010 0.015 

1/T (K' 1 ) 



Perdeutero 
0.020 0.025 



Fig. 5 Arrhenius plot of \nk(T) versus \/T comparing hydrogen 
and deuterium substituted donor- acceptor system. Solid line 
represents the contribution from off-diagonal electron-phonon 
couplings {gabn} °f a H the vibrational mode between eigenstate a 
and b which have frequencies in the range of 5 — 3200cm -1 , and 
dotted line represents the contribution from off-diagonal 
electron-phonon couplings of those vibrational modes which have 
the frequencies < 2000 cm~ l i.e. dotted line: g a t, n < 2000 cm~ l . 



in which our transformed Hamiltonian becomes 



3 Non-adiabatic relaxation rates 



Having parameterized our model, we can move on to compute 
the electronic energy transfer rate between the LPPP5 and PDI 
parts of the molecule. We now take the off-diagonal (non- 
adiabatic) electron-phonon coupling terms as the weak pertur- 
bation in order to describe transitions between the electronic 
eigenstates. In Ref. 20 Pereverzev and Bittner developed a 
time-convolutionless approach for treating non-adiabatic tran- 
sition for systems described by the Hamiltonian given above. 
We begin by performing the polaron transform using 



#0 = U~ l H U 

a 

+ I^^+to^ + W' 

where the renormalized electronic energies are 



t^a — t^a 7 . 



r 2 

K - r aan 



Applying the same unitary transformation to V gives 

V = Y,\a)(b\M abn , 



(14) 



(15) 



(16) 



(17) 



aim 



u = 



-Lanj$^\a)(a\(ai-a n ) 



Z\a)(a\e-^d-a») (U) 



where the system-bath operators are 

y ( t 2Gaan\ Y ^ G ^-G bbm ) ( J 



j . (18) 
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where 



1000 2000 

cj (cm -1 ) 



1000 2000 

cj (cm -1 ) 



(as 



(b) 



Fig. 6 Diagonal (a) and off-diagonal couplings (b) as in Fig [3] 
except for fully deuterated species 



At this point it is useful to connect the various terms in our 
Hamiltonian with specific physical parameters. The terms in- 
volving (Gaan — Gbbn)/h®n can be related to the reorganiza- 
tion energy 

where S n is the Huang-Rhys factor for mode n which is related 
to the Franck-Condon factor describing the overlap between 
the v n — 1 vibronic state in one electronic state with the v n = 
vibronic state in the other. Likewise, the energy difference be- 
tween the renormalized energies is related to the driving force 
of the transition, 



AE, 



ab '■ 



-E h . 



In the transformed picture the electronic transitions from state 
I a) —> \b) are accompanied not only by the creation or annihi- 
lation of a single phonon of mode n but also by the displace- 
ments of all the normal modes. Once we have transformed 
to the shifted (polaron) representation, we can construct the 
golden-rule rate in terms of the autocorrelation function 



K ab 



■L 



2Re / dxC ab (x)e 



-i(e b -e a )t 



(19) 



where k a b is the non-adiabatic rate (£jva)> and C a b(x) = 
(M a b(0)Mba(t))th is the autocorrelation of the polaron trans- 
formed electron-phonon operators 



Cab{*) 



Y^Gabn (Aabn(n n + l)e m " z - n n e- i<a ^) + O, 



'abnj 



n J 

X qabWfabi*), 



^■abn — 
^^abn 



Gbbn 



h(ti n 
(Gaan + Gbbn) 



and 



h®n 



1 



exp(h® n /k B T) - 1 ' 
is the Bose-Einstein occupation of the nth phonon mode. 



(21) 

(22) 

(23) 



Note that in the first term in Eq. 20 gives the non-Condon 
contributions to the correlation function. The two time depen- 
dent quantities are given by 



Gobi?) = e 



iLn( 



h&n 



sin(co n T) 



and 



Jab 



(x) = ^-2I,(^ + i)(^^) 2 (l-cos(co,T)) > 



(24) 



(25) 



These are related to the spectral density of the diagonal terms 
and as such are related to the Franck-Condon overlap integrals 
between the phonon modes in the a and b states. The decay of 
(M a b(0)Mb a (x)) t h is related to the decoherence time between 
the two electronic states and we have found that for the sys- 
tem at hand, the golden rule limit is clearly reached in 30 fs 
for the lowest temperatures considered in this study. Because 
our formalism does not include dissipative terms, there is a fi- 
nite recursion time for the correlation function. However, for 
larger systems this recursion time is very long and we assume 
that all correlation functions vanish for all time once the initial 
decay has occurred. 

In Fig. [5] we show an Arrhenius plot of the non-adiabatic re- 
laxation rates (k NA ) between the LPPP5 and PDI. At T =150K, 
kNA = 4.4ps _1 which is consistent with energy transfer rates 
observed in related systems. As T — » OK, kNA — > 1.6ps _1 re- 
flecting the contribution of zero-point vibrational motion. At 
higher temperature, we see log^A <* 1/T and conclude that 
the energy transfer is activated process with Ea = 7.56meV. 
In order to assess the relative contribution to the non-adiabatic 
rate from the two types of coupling modes identified in Fig. [3] 
we performed a series of rate constant calculations in which 
' we limited the non-adiabatic couplings to include only modes 
above or below a 2000 cm -1 cut-off. This cut-off lies clearly 
between the two peaks seen in the coupling spectra. Over the 
entire thermal range, there is no quantitive difference between 
the fully coupled model and a model including only couplings 
(20)from the C-H wagging modes. 
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Table 2 Kinetic data for energy transfer in the high-temperature limit (E a ). 



LPPP5-PDI-/* 



LPPP5-PDI-J 



coupling 



E A (meV) k(150K) ps" 1 E A (meV) k(150K) ps" 



§abn 



7.69 



4.19 



15.92 



1.74 



2 X gabt 



7.69 



16.77 



15.92 



6.99 



0.5 x g abi 



7.69 



1.04 



15.92 



0.43 



gabt cut-off below 2000 cm 



7.68 



4.19 



14.61 



1.72 



3.1 Kinetic isotope effect 

Having identified that C-H wagging modes drive the non- 
adiabatic energy transfer from the LPPP5 to PDI parts of the 
molecule, k^A should be sensitive to isotopic substitution of 
H for D. There are two factors contributing to an isotope ef- 
fect. First, the coupling strengths themselves depend upon 
the phonon frequency through G a bn = g a bn\/^/^^n- Since 
the g a bn do not depend upon the phonon frequency, we can 
roughly estimate that G ah ' n « 2 _1//4 G^ for modes involving 
proton motions. Since the correlation function and hence the 
golden-rule rate goes as G 2 abn , would give a factor of ~ y/l to 
the overall rate. The other contributions to the rate stem from 
the Franck-Condon terms arising from the displacement be- 
tween the two adiabatic potentials. Upon isotopic substitution, 
there will be shift in the zero-point energies of each state and 
generally, the nuclear wave function will be more localized. 
For the case at hand, the normal modes involving the proton 
motions give nearly all the contributions to the off-diagonal 
coupling (Gabn) but contribute little to the Franck-Condon fac- 
tors. Consequently, we can make a "back of the envelope" es- 
timate that ky 1 ) /k^ is approximately \/2 for systems where 
non-adiabatic transitions are mediated by proton motions. 

To test this, we determined the normal modes for 
perdeutero-LPPP5-PDI system in which all 35 protons were 
replaced by deuterons. Since the electronic transition moment 
and nuclear geometries for the ground and excited states do 
not depend upon the nuclear mass, we can use this data to 
construct a non-adiabatic model for energy transfer for the 
deuterated system. The coupling spectra for the perdeuterated 
system are shown in Figj6j Comparing to the perprotonated 
system, the g aan diagonal coupling terms are virtually iden- 
tical. This is not surprising since the g aan and g b bn terms in 
this system reflect C=C bond distortions and do not involve 
the C-H modes. This suggests that the Franck-Condon factors 
between the donor and acceptor states are more or less identi- 
cal for the perproto- and perdeutero- molecules. On the other 
hand for the off-diagonal non-adiabatic terms, g a b n , there is 
a systematic shift towards lower frequencies of all the major 
peaks in Fig. [6^ relative to those in Fig. |3]reflecting the effect 



of isotopic mass on the C-H frequencies. (co# « V^G)/)). 

In Fig. [5] we show the kinetic data for the fully deuterated 
LPPP5-PDI-d 35 and in Fig. we plot JfcW/ifcW over the 25K 
to 400K temperature range. Over the entire temperature range 
considered here, k^ /k^ > \/l which is consistent with our 
estimate above. In the high-temperature limit (T > 250^), 
we can estimate the Arrhenius activation energy as in the 
LPPP5-PDI-/*35 case. Comparing the Ea for the perproto- and 
perdeutero-cases, E^ ~ 2E\ \ In Table ^we given a sum- 
mary of the activation energies for various cases. Doubling or 
reducing the off-diagonal couplings by a factor of 2 has little 
effect on the activation energy indicating that that even if there 
were inaccuracies in the couplings themselves, the overall en- 
ergy transfer kinetics is determined by the proton vibrational 
frequencies. Finally, we considered the effect of truncating the 
coupling spectrum to include only the C-D wagging modes. 
As shown in Fig.|5j there is some sensitivity to the overall rate 
at higher temperature, but overall the kinetics are insensitive 
to this cut-off indicating that the C-D modes are the coupling 
modes in this case. 



4 Discussion 

In this paper we have presented a methodology for deter- 
mining non-adiabatic electronic transitions in large multi- 
chromophore systems starting from a common quantum 
chemical description. We believe this approach will be partic- 
ularly useful for studying systems in which one can not make a 
clear separation between the donor and acceptor species and in 
cases where the are common sets of vibrational motions. One 
crucial theoretical observation that comes out of this study is 
that in systems such as this, nuclear motions that are involved 
in the geometric reorganization of the molecule in a given ex- 
cited state may not necessarily be the same set of motions 
that are involved in the coupling between excited states. In 
fact for the system at hand, we identified that the C-H wag- 
ging motions couple the transition between the D and A state 
and not the C=C modes or ring torsional modes, even though 
these modes are involved in the geometric relaxation of the 
two states. We can understand this coupling in the following 
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Fig. 7 Ratio of non-adiabatic transition rates (k( H ' /k^ ) over the 
entire temperature range considered herein. 



way: In the absence of phonons, the two states are coupled 
via a transition moment, ju a b. Consequently, a photon with 
polarization along ii a j, induces an oscillation of the electronic 
state (of which the transition density would provide a snap- 
shot). In the same sense, in an IR transition the transition in- 
tensity is determined by the extent to which the permanent 
dipole moment of the molecule is modulated by a given nor- 
mal mode. Consequently, those modes with the largest pro- 
jection along the electronic transition moment will provide an 
oscillating electric field polarized in the right way to promote 
the transition. This is analogous to the situation in resonance 
Raman spectroscopy where vibrational modes along the tran- 
sition dipole are also coupled to the photon field and can par- 
ticipate in the electronic transition. Our identification of the 
nonadiabatic coupling vector with the state to state transition 
dipole moment is best considered a "propensity" rather than a 
selection rule. None the less it suggests that modes that are in 
the same irreducible representation as the dipole operator will 
give the strongest contribution to the non-adiabatic coupling. 

We also predict that in system such as this that the should 
be a clear kinetic isotope effect over the entire temperature 
range. While kinetic isotope effects are well known to occur 
in H transfer reactions 21 . In this case, changing the mass of 
the reacting atom modifies the zero-point energy and hence 
increases the activation energy required for bond cleavage. 
Here, we propose that the non-adiabatic energy transfer rates 
will be sensitive to nuclear mass. This effect is most anal- 



ogous to non-adiabatic relaxation dynamics of the solvated 
electron. In this case, ultrafast pump-probe experiments by 
Barbara's group indicate a k^/k^ = 1.4 in comparing the 
ground-state recovery of an excess electron in water versus 
heavy water^J. The effect was largely attributed to the quan- 
tum mechanical librational modes of the surrounding solvent 
cavit y 23 " 27 i 

The close analogy between non-radiative transition and en- 
ergy transfer events came be traced to early works by Robin- 
son and co-workers 28 " 32 ! and kinetic isotope effects in non- 
radiative relaxation rates for various deuterated aromatic sys- 
tems has been known since the 60 , sl 3Q l 32 l For example, deuter- 
ation dramatically increases the triplet lifetime of various 
deuteronapthalene systems as observed in ESR experiments 
by Hutchenson 33 . Similar effects have been reported for the 
3 B\ U — » l A\ g triplet to singlet conversion in benzene and 
benzene-^6 with kS H > /k^ D ' ranging between 2 and 4^1 For 
example, in Ar matrix at 4K, k^/k^ = 5.9EE3 while in 
EPA glass at 77K JfcW/ifcW = 1. 1 34 . It should be pointed 
out that theoretical studies of the isotope effect in the non- 
radiative decay of benzene have shown little agreement with 
experimental results. For example, in Ref. — , Burland and 
Robinson use a level counting approach and obtain an estimate 
of k^/k^ « 3.7 x 10 5 . Other theoretical work using saddle 
point methods to numerically evaluate the rate constant in the 
statistical limit give ratios of k^/k^ = 24 36 to k^/k^ = 
1.1 x 10 4 and 2.2 x 10 2 depending upon the parameters for 
the e2 g mode 37 . Finally, recent work by Zamstein using an 
ab initio parameterization and a phase-space approach for the 
non-radiative 3 B\ U — )■ l A\ g triplet to singlet conversion in ben- 
zene and benzene-^ gives a k^/k^ = 3.7 — 7.8^51. 
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